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INTRODUCTION 


V 

The  purpose  of  this  paper  is  to  survey  the  one-sided  search  prob¬ 
lem,  starting  with  Koopman's  work  in  1946  and  continuing  to  the 
present.  The  paper  summarizes  some  basic  results  for  both  the 
optimal  allocation  of  effort  problem  and  the  best  track  problem  for 
stationary  and  moving  targets;  but,  it  is  intended  as  an  introduc¬ 
tion  to  the  literature  of  the  field  rather  than  an  exposition  of 
it. 


The  paper  does  not  include  two-sided  searches,  l.e.,  searches  in 
which  the  target  reacts  in  any  intelligent  way  to  the  searcher.  It 
does  not  include  surveillance  problems  or  problems  involving  false 
contacts  or  decoys.  It  includes  only  repeatable  searches. 


This  survey  was  aided  by  earlier  surveys,  particularly  J.  M. 
Dobbie'a  published  in  Operations  Research  in  1968,  and  Marc 
Mangle's  1980  OEG  publication.  The  bibliography  includes  not  only 
those  papers  actually  referenced  in  the  present  work,  but  also 
other  papers  of  interest  in  the  area  of  search. 


THE  SEARCH  PROBLEM 


Search  theory  is  one  of  the  oldest  areas  of  operations  research. 
Problems  involving  search  arise  in  such  diverse  areas  as  the 
military  looking  for  enemy  submarines,  the  Coast  Guard  searching 
for  small  boats  lost  in  a  storm,  prospectors  surveying  for  mineral 
deposits,  the  forest  service  looking  for  missing  backpackers,  law 
enforcement  officers  searching  for  lost  weapons  or  escaped  crimi¬ 
nals,  a  secretary  looking  for  a  missing  file,  or  an  analyst  scan¬ 
ning  a  computer  printout  for  a  particular  piece  of  data.  All  of 
these  problems  have  two  elements  in  common  --  a  target,  in  the 
broad  sense  of  something  being  searched  for,  and  a  searcher. 

There  are  usually  two  types  of  cost  involved  in  search  problems. 

The  first,  the  cost  of  the  search  itself,  may  be  measured  in  such 
terms  as  dollars,  time,  manpower  expended,  or  fuel  expended.  We 
often  want  to  search  in  such  a  way  as  to  maximize  the  probability 
of  finding  the  target  at  a  minimum  cost,  or  until  our  resources  run 
out  (fixed  cost).  A  second  cost  is  the  cost  of  not  finding  the 
target.  This  cost  may  be  measured  in  dollars,  in  inconvenience,  or 
even  in  lives  lost.  The  two  costs  need  to  be  balanced  in  each 
search  situation.  In  general,  we  want  to  devise  a  search  plan,  or 
"track"  which  uses  the  resources  most  effectively  under  such 
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physical  limitations  as  the  terrain,  the  searchers,  the  instruments 
used,  the  resources  available,  and  the  nature  of  the  target  itself. 

The  search  problem  can  be  loosely  described  sb  follows: 

•  The  target  is  located  in  an  area  which  Is  much  too  large 
for  the  searchers  to  search  completely. 

•  The  location  of  the  target  is  not  known  exactly,  but 
probabilities  can  be  associated  with  subregions  of  the  main 
search  area. 

•  The  target  may  (or  may  not)  move. 

•  One  or  more  searchers  may  look  for  the  target,  and  they  may 
use  detection  equipment  to  do  so. 

In  order  to  solve  the  problem,  we  need: 

•  A  model  for  the  location  of  the  target  at  the  start  of  the 
search.  This  model  we  will  call  the  initial  density. 

•  A  model  for  the  motion  of  the  target,  which  we  denote  by 
q(x,y, t). 


3 


•  A  mathematical  goal,  or  objective  function,  such  as  mini¬ 
mizing  the  time  to  find  the  target,  or  maximizing  the  prob¬ 
ability  of  detection  by  time  T, 

The  search  literature  breaks  the  problem  down  into  two  main  cate¬ 
gories:  optimal  allocation  of  effort  problems,  and  best  track 

problems*  The  optimal  allocation  of  effort  problem  is  mathematic¬ 
ally  the  easier  of  the  two  problems  and  hence  more  work  was  done  in 
this  area  earlier.  Optimal  allocation  of  search  effort  may  mean, 
for  example,  optimizing  the  amount  of  time  spent  searching  in  each 
subarea.  These  problems  are  nice  mathematically  since  it  is  often 
possible  to  prove  that  a  plan  is  optimal.  However,  it  may  be  that 
the  optimal  plan  is  not  "doable."  For  example,  in  the  figure 
below,  the  plan  may  say  to  put  50  percent  of  the  total  effort  in 
(3,2)  and  50  percent  in  (3,4).  Due  to  physical  constraints,  this 
allocation  may  be  impossible, 
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e.g.,  the  searcher  may  not  be  able  to  get  from  (3,2) 
without  expending  some  effort  in  (3,3).  It  would  be 
for  the  searcher  to  have  a  track  to  follow.  But  the 


to  (3,4) 
more  useful 
best  track 
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problems  are  difficult  to  solve*  Work  has  progressed  significantly 
only  since  the  mid-seventies. 

We  will  look  first  at  the  best  allocation  of  effort  problem  and  its 
history. 


KOOPMAN ,  1946 


The  earliest  developments  in  search  were  made  by  Bernard  Koopman 
and  his  colleagues  In  the  Anti-Submarine  Warfare  Operationa 
Research  Group  (which  later  became  the  Operations  Evaluation  Group) 
of  the  U.S.  Navy  during  World  War  II.  Their  purpose  was  to  aid  the 
Navy  in  finding  efficient  ways  to  search  for  enemy  submarines.  The 
work  done  from  1942  to  1945  was  published  in  a  book,  Search  and 
Screening  (Koopman,  1946).  Originally  classified  Confidential,  the 
work  was  declassified  in  1958.  A  new  edition  of  Search  and 
Screening  was  published  in  1980.  Many  of  the  results  from  the  OEG 
work  were  published  in  a  aeries  of  articles  in  Operations  Research 
in  the  mid-fifties  (Koopman,  1956  a&b,  1957).  Koopman's  work  is 
basic  to  search  theory  and  is  a  good  place  to  start  a  survey  of 
this  field. 

Let  us  look  first  at  the  law  of  random  search*.  Suppose  there  is 
a  region,  A,  over  which  a  search  must  be  made.  Assume  nothing  is 
known  about  the  location  of  the  target  except  that  it  ia  in  A, 
i.e.,  we  will  assume  a  uniform  target  distribution  and  that  the 
target  is  stationary  relative  to  the  searcher.  Assume  that  if  we 
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See  Koopman  (1946),  p.28  or  Koopman  (1980),  p.  71 


pass  within  W/2  distance  units  from  the  target,  we  will  detect  it 
with  probability  1.  Assume  also  that  the  searcher  takes  a  random 
piece-wise  linear  path  of  total  length  L,  and  l*  L/n  ia  the  length 
of  one  of  n  equal,  rather  long  (in  relation  to  W)  segments. 


Koopman  shows  that  the  probability  of  finding  the  target  under 
these  conditions  is 

P  .  1  - 

Since  the  search  path  was  random,  thiB  probability  is  usually  lower 
than  it  would  be  if  we  had  searched  systematically.  P  is  then  a 
lower  bound  for  the  actual  probability  of  detection  and  is  useful 
for  approximation  purposes.  If  we  assume  that  it  is  only  necessary 
to  search  in  a  straight  line  path,  P  becomes 

P  -  WL/A. 
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This  estimate  of  P  assumes  "ideal"  conditions  and  thus  provides 
an  upper  bound  for  the  probability,  as  indicated  in  the  figure 
below , 


The  true  probability  lies  somewhere  in  the  hatched  region. 

Koopman  looked  at  the  optimal  allocation  of  effort  problem  in  both 
discrete  and  continuous  space.  HIb  objective  was  to  optimize  alio 
cation  of  search  effort.  Suppose  that  we  have  a  search  area 
divided  into  two  parts  A^  and  A2*.  Suppose  that  the  prob¬ 
ability  that  the  target  is  in  A^  iB  p^,  i  ■  1,2,  and  define 
,  the  density  of  search  effort  in  A^  by 

Effort  in  A, 


*See  Koopman  (1946)  p.  35,  or  Koopman  (1980)  p.  146. 


Koopman  found  that,  assuming  the  target  is  stationary  (or  more 
generally,  stationary  relative  to  the  searcher)  and  assuming  the 
lav  of  random  search,  one  of  three  situations  arises: 

If  Pj/Ai  is  much  larger  than  P2/A2,  aearc^  only  A^. 

If  P2^2  *s  muc*1  larger  than  P^/A^,  search  only  A2> 

Otherwise,  distribute  the  total  effort  A^^  +  A2d>2 

using  4>l  "  log(P1/Ai)  -  (1/A)  [Aj^logfPj/Aj)  +  A  2log  (P  2/A2 )  ]  +  4>/A. 

For  the  continuous  case,  Koopman  defined  the  problem  as  follows: 

•  The  target  is  stationary. 

•  It  is  contained  in  a  region  A.  The  probability  before  the 
search  starts  that  the  target  is  in  (x,x+dx)  and  (y,y+dy)  is 
p(x, y )dxdy . 

«  The  total  search  effort  is  <t>. 

4>(  )  is  the  density  of  search  as  a  function  of  x.  If  we  integrate 
<t>(.  )  over  any  subregion  B  of  A,  we  get  the  effort  expended  in  B. 
Koopman  shows  that  the  problem  reduces  to  maximizing 
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p<*)  -// ptx.yXl  -  e'^  (x'y))dxdy 
A 

t 

where  p(x,y)  Is  the  probability  that  the  target  la  around  (x,y) 

and  1  -  e”*<x’y)  la  the  probability  of  finding  it  there  assuming 
random  search.  Using  Lagrange  multiplier  and  Calculus  of  Varia¬ 
tions  techniques,  he  finds  the  unique  b  such  that 

log(b)  -  C 1 / )//log  p(x,y)  dxdy  +  <t>/Ab  ■  0 
Ab 

where  A.  is  a  subset  of  A  which  depends  on  b.  Then  $(x,y)  " 

D 

loglp(x, y)/b]  is  the  function  defined  by  //*(  x,  y  )dxdy  ■  WL  -  4> 

A 

and  (x,y)  >.  0  which  gives  P(^)  its  largest  valu<£.  Geometrically, 
we  can  find  this  value  of  b  by  first  plotting  z  -  log  p(x,y) 
over  the  region  A.  Cut  this  surface  by  the  horizontal  plane 
z  ■  log  b  in  Buch  a  way  that  the  volume  of  the  surfaces  above  the 
plane  is  equal  to<J>.  The  orthogonal  projection  on  the  x,  y  plane 
will  be  A  . 

D 

For  both  the  discrete  and  continuous  cases,  Koopman  observed  that 
for  this  problem  effort  allocation  has  an  additive  property.  That 
is,  if  <f>  units  of  effort  are  available  at  the  beginning  of  the 
search  and  4>*more  units  become  available  later,  the  best  search 
plan  remains  the  one  which  would  have  been  chosen  at  the  beginning 
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had  we  known  the  total  effort  would  be  $  +  «  In  other  words, 

says  Koopman 


"A  well  planned  search  cannot  be  improved  by  a 
redistribution  of  search  made  at  an  intermediate 
stage  of  the  operation  in  an  attempt  to  make  use 
of  the  fact  that  up  to  that  time  the  target  had 
not  yet  been  observed.'*1 


koopman  (  1  946),  p.38  or  Koopman  (1  980),  p .  151. 
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SUBSEQUENT  DEVELOPMENTS 


Variations  on  both  of  Koopmart's  problems  followed  quickly  after  the 
publication  of  the  unclassified  version  of  his  work.  The  discrete 
case  easily  extends  to  n  boxes  and  to  many  variations. 

Blachman  (1959)  looked  at  the  problem  of  finding  an  object  in  one 
of  n  boxes.  He  assumed  that  the  probability  the  object  appears 
In  each  box  is  known,  and  the  time  of  appearance  is  uniformly 
distributed  over  a  large  interval.  Blachman  and  Proschan  (1959) 
derive  an  optimum  search  procedure  for  a  generalization  of 
Blachman's  problem.  In  this  case,  the  object's  arrival  1b  in 

accordance  with  a  Poisson  process  with  arrival  rate  A.  A  cost  of 

looking  in  each  box  Is  added.  They  consider  which  boxes  to  scan 
and  how  to  schedule  the  scans  to  minimize  the  time  between  the 
arrival  of  the  object,  and  finding  it.  The  model  Is  applicable  to, 
say,  the  appearance  of  miasiles  for  which  early  detection  is 
crucial. 

Gilbert  (1959)  added  to  Koopman's  two  cell  search  a  cost  (perhaps 
in  time)  for  switching  from  one  cell  to  another.  He  treats  the 
problem  as  a  kind  of  one  person  game,  and  applies  the  method  to  a 

search  for  an  odd  sized  bolt  In  one  of  two  bins.  Gluss  (1961) 

looked  at  the  r-box  problem  if  the  boxeB  are  all  in  a  .line  and  Ross 
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(1968)  added  a  reward  Rj  for  finding  the  target  in  box  i. 

Wegner  (1980)  find  necessary  and  sufficient  conditions  for  the 
existence  of  admissible  search  strategies  which  minimize  the 
expected  cost  of  "at  least  Burely"  finding  an  object  when  overlook 
probabilities  are  Included.  He  gives  r  procedure  for  computing  an 
optimal  strategy* 

Most  of  the  above  problems  were  solved  by  a  dynamic  programming 
approach.  Matula  (1964)  derived  conditions  for  the  existence  of  an 
ultimately  periodic  search  with  minimum  cost.  He  finds  a  closed 
form  solution  to  this  problem  rather  than  the  recursive  solution  of 
dynamic  programming.  Pollock  (1964)  introduced  a  Bayesian  approach 
to  the  optimal  allocation  problem.  Decisions  are  made  sequentially 
based  on  what  had  been  observed  until  that  time  to  minimize  the 
expected  cost  of  searching  and  making  wrong  decisions.  In  1970, 
Pollock  (1970)  considered  the  case  when  the  target  is  moving 
between  two  regions.  The  target  moves  in  a  Markovian  fashion  and 
with  known  parameters.  He  tries  to  find  the  expected  number  of 
"looks"  required  to  find  the  target,  and  solves  the  problem  in 
certain  special  cases  by  means  of  dynamic  programming. 

The  preceeding  papers  all  deal  with  the  problem  of  effort  alloca~ 
tion  when  the  search  Bpace  is  discrete.  The  following  papers  use 
continuous  search  space,  or  both  continuous  and  discrete  space 
through  the  use  of  the  Stieltjes  integral.  Most  of  the  recent 
papers  use  this  latter  approach. 
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Charnes  and  Cooper  (1958)  showed  that  mathematical  programming 
could  be  used  in  search  problems  to  look  at  broader  classes  of 
problems.  They  applied  convex-programming  along  with  the  Kuhn- 
Tucket  conditions  to  obtain  the  solution.  The  algorithm  they 
obtained  was  an  important  step  in  solving  search  problems  by 
c  o  mp  u  t  e  r . 

In  1961,  deGuenin  generalised  Koopman's  models  in  an  algorithm 
which  made  no  assumptions  on  the  detection  probability  function. 

He  felt  that  the  law  of  random  search  (negative  exponential)  is  not 
valid  for  many  non-military  applications.  He  used  instead  func¬ 
tions  he  terms  as  "regular"  --  that  is,  the  graph  ia  strictly  con¬ 
cave  downward,  passes  through  the  origin  with  a  tangent  of  positive 
slope,  and  increases  monotonlcally  to  a  horizontal  aBympotote  no 
higher  than  positive  one.  Many  detection  functions  UBed  by  Koopman 
have  this  property,  and  most  subsequent  mathematical  developments 
have  been  based  on  deGuenin's  regular  functions. 

Dobbie  (1963)  started  with  the  additive  property  mentioned  by 
Koopman  and  used  this  property  to  derive  the  optimal  search  distri¬ 
butions.  He  was  the  first  to  define  passive  observations  --  the 
target  does  not  react  to  the  search,  the  search  does  not  materially 
change  the  target  (as  covering  a  lost  item  with  dirt)  or  the 
searcher  (as  with  fatigue)  --  as  opposed  to  active  observations  and 
to  discuss  the  mathematical  consequences.  He  pointed  out  that 
often,  especially  with  active  observations,  the  detection  function 
is  not  regular  (or  deGuenin). 
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Zahl  (1963)  derived  necessary  and  sufficient  conditions  for  the 
exiatance  of  solutions  to  the  problem  of  maximizing  the  detection 
probability  with  a  given  effort. 

Various  sequential  formulations  of  the  search  problem  have  been 
applied  to  several  fields.  Engel  (1957)  looked  at  the  search  for 
certain  minerals  as  a  two  stage  process.  The  first,  and  least 
expensive,  consists  of  one  or  more  preliminary  searches.  The 
second  search  is  detailed  and  more  expensive,  and  occurs  only  where 
the  required  "clusters"  of  signals  in  the  preliminary  search  indi¬ 
cate  the  target  is  likely  to  be  found.  Posnor  (1963)  uses  a 
similar  preliminary  scanning  technique  to  search  for  a  loBt  satel¬ 
lite  by  radar.  DeGuenin  (1963)  adds  a  middle  stage,  screening  the 
data,  when  searching  for  oil  wells. 

Stone  made  use  of  Calculus  of  Variations,  convexity  properties,  and 
generalized  Lagrange  multiplier  techniques  (which  allow  for 
inequality  constraints  and  do  not  require  differentiability  assump¬ 
tions)  to  Eormulate  a  systematic  treatment  of  Beurch  theory  in  his 
1975  book  Theory  of  Optimal  Search.  In  this  book  he  deals 
primarily  with  stationary  targets  but  extends  hit)  methods  to  false 
targets  and  Markovian  motion. 

Conditionally  deterministic  target  motion  motion  in  which  the 
initial  position  and  speed  of  the  target  is  known  but  Its  direction 
is  not  --  was  considered  by  Stone  (1973)  and  Pursiheimo  (1977). 


Mangle  (1980)  Included  an  algorithm  for  Markovian  target  motion  in 
which  the  moving  target  motion  is  reduced  to  a  sequence  of 
stationary  target  problems.  The  observation  that  a  search  plan 
maximizes  the.  overall  probability  of  detecting  a  moving  target  if 
and  only  if  it  maximizes  the  probability  of  detecting  a  stationary 
target  at  discrete  time  intervals  was  made  by  Brown  (1980).  When 
there  is  random  target  motion,  the  number  of  possible  target  paths 
is  infinite.  Stone  (1979)  worked  on  this  problem.  When  the 
detection  function  Is  concave,  he  gave  conditions  for  optimal 
search  plans  which  include  any  "reasonable"  target  motion. 

Stone,  et.  al.  (1978)  summarizes  the  optimal  allocation  problem  and 
included  several  algorithms  for  its  solution.  They  found  necessary 
and  sufficient  conditions  for  optimal  search  for  a  moving  target 
when  time  is  discrete  and  an  exponential  detection  function  is 
assumed. 

For  a  search  for  a  moving  target  to  be  optimal,  it  is 
necessary  and  sufficient  that  at  each  time  t  it  assigns 
an  allocation  which  1b  optimal  for  the  stationary  target 
problem  which  one  obtains  at  time  t  by  conditioning  on 
failure  to  detection  ufter  t  as  well  an  before  t  under 
the  plan.1 

The  algorithms  first  find  an  optimal  allocation  of  effort  for  the 
initial  target  distribution.  Then  for  times  t  »  2,...,T,  they 

one ,  et. .  al.  (1978),  abstract. 


16 


calculate  the  posterior  distribution  for  the  target  location  at 
time  t  given  failure  to  detect  at  all  previous  times  and  allocate 
effort  for  that  stationary  target  problem.  The  plan  resulting  from 
the  first  pass  is  called  the  "myopic"  plan.  Subsequent  passes  are 
made  reallocating  the  effort  each  time.  One  can  come  as  close  to 
the  optimal  plan  as  desired  by  performing  enough  passes.  In  many 
situations,  the  myopic  plans  are  almost  as  good  as  the  optimal 
plan. 


The  optimal  plain  typically  pay  a  penalty  in  probability 
of  detection  at  early  hours  in  order  to  maximize  that 

probability  at  time  t.  ...when  the  myopic  plan  is  close 
to  optimal,  the  myopic  plan  is  a  good  one  for  operational 
purposes. * 


*  St  one,  et.  al.  (1978),  abstract. 


KOOPMAN  REVISITED 


In  1979,  Koopman  (1979,  a&b)  published  two  papers,  and  in  1980 
republished  his  1946  work,  Search  and  Screening.  In  1979a  he 
generalized  his  1946  work.  In  1979b  his  emphasis  was  on  opera* 
tlonal  feasibility;  this  subject  will  be  discussed  in  the  last 
section  of  this  paper.  Both  of  these  papers  ate  incorporated  in 
the  1980  edition  of  Search  and  Screening. 

Koopman  rederived  his  law  of  random  search  under  Blightly  more 
general  assumptions.  He  assumed: 

e  Search  is  by  passive  observations. 

e  The  search  is  repeatable. 

e  Short  range  detectors  are  used. 

e  All  factors  (target,  searcher,  and  environment)  are 

constant. 

•  d>(x,y)  is  the  search  effort  applied  at  (x,y). 
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Define  D(x,ytO,  where  z  *  d>(x,y),  as  the  probability  the 
target  la  detected  given  that  it  ia  at  (x,y)  and  with  search  effort 
z.  Then 

D(x,y,z)  -  1  -  e  W<x,y)z,  w(x,y)  >  0 

In  general,  the  larger  w(x,y),  the  greater  the  probability  of 
detection}  it  ia  a  sort  of  local  measure  of  detectability  at  (x,y). 
Under  the  assumption  that  everything  within  a  range  R  of  the 
searcher  is  detected,  w  becomes  the  aweepwidth  defined  in 
Koopman's  1946  paper. 

The  optimal  distribution  of  searching  effort  problem  can  be 
formulated  as  follows. 

Find  the  search  density  function  $(x,y)  which  maximizes 

F(d>)  ■  ffp  (  x,  y  )D(x,  y  ,  d*(x,  y  ) )  dxdy 
A 

■  ffp(,x,y)  |l  “  exp  [-w(x,  y  )$(x,  y )  ]  j  dxdy 
A 

subject  to  //««  ,y)dxdy  ■  ,  and  $(x,y)  >_  0. 

A 

<t>is  the  total  available  effort. 
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He  assumed  further  that  the  functions  p,  w,  and  <^are  all 
continuous.  He  used  a  method  developed  by  Gibbs  in  1928  to  solve 
this  Calculus  of  Variations  problem.  Koopman  shows  that  there  is  a 
unique  X  (similar  to  the  "b"  in  his  earlier  work)  Buch  that  the 
optimal  ^ (x, y )  is 

$(x,y)  -  [l/w(x,y)I  log  (p(x»y  )w(x,  y)/  X  )  for  (x,y)«A(4), 
and 

*■//  I  log[ p (x, y  )v (x, y )  ]  -  log  X  |  dxdy/w (x, y >. 

A  (*) 

The  geometric  interpretation  is  the  same  as  in  the  earlier 
f or  mula  t ion. 

By  replacing  the  exponential  detection  function  D(x,y,z)  with  a 
partial  detection  function  b(x,y,z),  where 

00 

b(x,y,z)  »  1  -fe  WzdG(w) 
o 

(the  Integral  1b  a  Laplace-St ielt jes  transform),  Koopman  findo  a 
larger  class  than  the  regular  (deGuenln)  detection  functions  for 
which  his  results  are  valid.  G(w)  usually  involves  (x,y);  as 
w  goeB  from  0  to  po  ,  G(w)  goes  non-decrsasingly  from  0  to  l. 

(For  other  properties  of  G  see  Koopman  (1979a)).  However, 

Koopman  points  out  that  "The  practical  problem  of  finding  G(w) 
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has  so  far  only  been  solved  by  guesswork,  without  subsequent 
ver if i ca t ion.  "  *  The  Normal  distribution  has  been  suggested. 
Stone  (1975),  and  Richardson  and  Belkin  (1972),  have  studied  the 
use  of  the  gamma  distribution  when  the  sweep  width  is  uncertain. 


^Koopman  (1979b),  p.538. 
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BEST  TRACK  PROBLEMS 


V 

V( 

Optimal  effort  allocation  problems  are  relatively  easy  to  solve; 
however,  it  is  common  that  the  goal  ik'  to  find  the  best  search 
track.  Mathematically  this  is  a  harder  problem,  and  often 

approximate  solutions  are  all  that  a,re  available. 

V 

In  1946,  Koopman  developed  the  method  of  parallel  sweeps  for 
searching  for  a  stationary  target  or  a  target  whose  speed  and 
direction  are  known.  In  either  case  it  is  assumed  that  the  prior 
target  distribution  is  uniform  in  the  search  area.  This  search 
plan  calls  for  the  searchers  (e.g.,  airplanes)  to  move  along  a 
series  of  parallel  lines  (whose  distance  apart  depends  on  the 
search  environment)  which  cover  the  area.  In  the  case  of  a  moving 
target,  the  sweeps  are  made  parallel  to  the  target's  motion. 

Between  the  1940s  and  1970  almost  no  progress  was  made  in  solving 
best  track  problems.  In  1974,  Lukka  (1974)  worked  out  the  theory 
of  optimal  track  for  stationary  targets,  targets  whose  motion  is 
known,  and  targets  whose  motion  is  almost  known.  The  methods  rely 
on  the  theory  of  optimal  control.  Mangel  and  Thomas  (1979)  wrote  a 
tutorial  type  paper  developing  from  first  principals  analytical 
methods  for  search  for  a  moving  target. 


22 


Mangel  (1981),  basing  his  work  on  Lukka's,  derived  algorithms  for 
the  optimal  control  result,  one  where  the  detection  rate  is  inde¬ 
pendent  of  velocity  and  one  where  it  is  not.  Mangel  defines 
f(x,t,z)  and  u(x,t,z)  as 

f(x,t,z)dx  "  Prob[x(t)  <  (x,x+dx)  and  search  along 
z(r),  0  <  t  was  not  successful] 

u(x,t,z)  ■  Problnon-detection  to  time  T|x(t)  ■  x, 
and  search  along  z(r),  t  <r^  T] 

He  showed  (Mangel  and  Thomas,  1979)  that  these  quantities  must 
satisfy 

3f  -  s  1/2  a2  (a  f)  -  £  d  (b.f)  -  ^  (x.  t,z)f  m 

3tT  i,j  j  ■  i  9x1 

with  f(x,0,z)  ■  fo(x)  plus  boundary  conditions, 

du  -  -  2  1  /  2a .  ,  ~  £  b,  Su  +  *  (x, t,z)u  (2) 

ij  i  ^7 

with  u(x,T,z)  ■  1  +  boundary  conditions, 

where  is  the  diffusion  matrix  and  b^(x)  represents  the 

average  velocity  of  the  target  (the  drift  coefficient).  \^(x,t,z) 
is  the  instantaneous  detection  rate,  defined  by 
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^(x,t,z)dt  •  Prob[detection  in  (t ,  t+d t  )|  X ( t  )  ■  x, 

Z(t)  ■  z,  no  detections  before  t] 

Mangel  modified  the  ray  techniques  of  J.B.  Keller  (Keller,  1978)  to 
solve  equations  (1)  and  (2). 

To  find  the  best  search  track  which  maximizes  P(t),  the 
probability  of  detection  at  time  T,  we  must  pick  v(t),  the 
searcher  velocity  so  that  P(T)  ■  1  -/f(x  ,T,z)dx  is  maximized, 
dz/dt  ■  v(t),  z(0)  m  zq,  and  £(x,T,z)  satisfies  (1).  Lukka 
(197  4)  derived  the  maximum  principle  for  this  problem,  showing  that 
the  optimal  velocity  v*(t)  makes  H(t,v,A)  a  maximum  over  all 
allowed  velocities,  where  f(x,t,z)  satisfies  (1),  u(x,t,z) 
satisfies  (2),  dz/dt  ■  v(t),  and 

H  (t  ,  v ,  A )  -  ff  (x,  t ,  z  )u  (x,  t ,  z  )^(x,  t  ,z  )dx  fZA^tJv^t), 

dA  /dt  -  ~f  f  (x,  t ,  z  )u(x,  t ,  z  )4^-(x,  t ,  z  )  dx,  with  A(T)  ■  0. 
i 

Mangel,  points  out  that  when  the  detection  function  is  independent 
of  velocity,  the  searcher  should  move  at  maximum  speed  and  in  the 
direction  of  A(X)  where 


A  ( t  )  -  A  (0  )  -//f(  x,  s,  z  )u  (x,  .>i,  z  )^dxds . 
0 
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"‘•nr  vn  i  — -  - 


A  (0  ) 


f(x,t,z)u(x,t,z )^dxdt. 


It  seems  clear  that  we  are  a  long  way  from  the  routine  solution  of 
problems  of  this  sort. 


OPERATIONAL  IMPLEMENTATION 


In  his  "Operational  Critique  of  Detection  Laws,"  Koopman  asks, 
"Given  a  theoretically  perfect  solution  to  a  problem  of  optimal 
search,  how  accurately  can  it  be  implemented  by  the  dispositions  of 
paths  of  real  searchers?"1  His  1946  work  looked  at  this  question 
for  several  situations  and  decided  that  between  70  and  85  percent 
of  the  theoretical  optimum  is  the  best  which  can  be  expected* 

Since  this  is  the  cose,  he  pointed  out  that  it  may  be  better  to 
spend  time  finding  "good  enough"  solutions,  i.e.,  finding  good 
"useful  first  approximations"  rather  than  continuing  to  find 
elegant  exact  solutions. 

Many  of  the  techniques  discussed  in  this  paper  require  extensive 
computer  time  to  implement.  A  solution  which  requires  an  hour  on  a 
larger  computer  is  not  very  useful  for  searchers  operating  from  an 
aircraft  carrier  or  from  a  small  law  enforcement  office.  The 
searchers  need  to  know  in  real  time  where  to  look  and  how  long  to 
look  there.  In  order  to  make  the  theory  more  accessible  to  those 
who  need  it,  much  work  is  needed  to  devise  good,  simple  approxima¬ 
tions  to  optimal  plans. 


koopman  (1979a),  p.131. 
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